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1. Introduction 

The most notorious and renowned probability distribution (PD) used in the field 
of statistical mechanics is undoubtedly that deduced by Gibbs for the canonical 
ensemble 012], usually referred to as the Boltzmann-Gibbs (BG) equilibrium dis- 
tribution 

(!) PgW = ? , 

^BG 

with Ei the energy of the microstate labeled by i, j3 = 1/fceT the inverse tem- 
perature (T), ks Boltzmann's constant, and Zbg the partition function. The 
exponential term Fbg — exp (—(3E) is, of course, called the BG factor. Recently 
Beck and Cohen E] have advanced a generalization (called "superstatistics" ) of 
this BG factor, assuming that the inverse temperature j3 is a stochastic variable. 
They effect a multiplicative convolution to obtain a generalized statistical factor 
Fqs hi the fashion 0] 



/>oo 

(2) F GS = / d(3 (3f(f3) exp {-(3E) = f o F, 

Jo 

where /(/3) satisfies 

/■oo 

(3) / d/?/G9) = 1. 



BG, 



Note that the o— sign is used to denote the multiplicative convolution between two 
PDs: multiplicative convolution of two PDs fx and fy is the PD of the product 
of the two corresponding random variables X and Y, Superstatistics would arise 
as a "statistics of statistics" that takes into account fluctuations temperature |3J 
0] El E| 03 • Fixed temperature values, as those emblematic of the canonical 
ensemble of Gibbs', require that the system to be described be in thermal contact, 

l 
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and in equilibrium, with an infinite heat reservoir characterized precisely by that 
temperature T. 

Beck and Cohen also show that, if f(j3) is a x 2 distribution, nonextensive thermo- 
statistics (NEXT) is obtained, which is of interest because NEXT, also known as 
Tsallis' thermostatistics, is today a very active field, perhaps a new paradigm for 
statistical mechanics, with applications to several scientific disciplines jl()l 1111 H^] . 
In working in a NEXT framework one has to deal with power-law distributions, 
which are certainly ubiquitous in physics (critical phenomena are just a conspicu- 
ous example Q2]). Indeed, it is well known that power-law distributions arise quite 
naturally in maximizing Tsallis' information measure (q is a real parameter called 
the "nonextensivity index" ) 

(4) H q (/) = ^—(l- | + °° f{xYdx 

subject to appropriate constraints. In the case of the canonical distribution there is 
only one constraint, the energy E, i.e., (E) = K (K a constant) and the equilibrium 
canonical distribution writes f(x) = (l/Z q )[l — (1 — q)f3 q E] i- 1 , with f3 q and Z q 
standing for the NEXT counterparts of (3 and Zbg above. It is a classical result 
that as q — > 1, Tsallis entropy reduces to Shannon entropy 

/+oo 
f(x)logf(x). 
-oo 

Tsallis' PDs are encountered in analyzing a rather large variety of physical sys- 
tems ^] El El El H8| . which encourages people to continue investigating the 
nonextensive formalism along multiple viewpoints and a multitude of paths. 

In such a spirit, and further pursuing along the road first travelled by Beck and 
Cohen P], we show here that a new type of superstatistics can be constructed for 
which the fluctuating physical quantity is the nonextensivity index q. 

2. Mathematical framework 
2.1. Notations. Using notation x+ = max(x,0), let us denote 

(6) S q ^ (x) = A q ( 1 - ^ 

the Tsallis distribution with variance a 2 ^ q \ , nonextensivity index q, and normal- 
ization constant 



(7, A.- F( - + 3/2 



r(^r + i) 

A word of caution is necessary here. There are several versions of Tsallis' thermo- 
statistics. In some of them Tsallis' PDs carry the power 

(i) while for others it is 

(8) (ii) (our choice above). 

Of course, translation between the two is straightforward. If we denote (a) by q* 
the nonextensivity index of the first choice above, and (b) by q that of the second 
one, then 



(9) 



q* = 2 - q. 
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2.2. Main result. Before stating our main result under the form of Theorem^ 
we need the following lemma. 

Lemma 1. Vc, Va > 0, there exists a discrete distribution pk (a) such that 

( 10 ) _ _ 



J^Pk (a) S 1+ i 2aa 2 (a 



k=0 



crf(y / a)c 





else 



Proof. Assuming x G [— er\/2a, +a^/2a\ , let us compute the left hand side of I|1U|) 
as 



+oo 



fc=0 



a v 2a 



r(fc + 3/2) 



exp(-a) ^ a fc+1 / 2 ,. _ , { 

erf(v^) ^ T(/c + 3/2) a V2aT(k + IJ^tF V 2 a( r 2 

exp(— a) 1 

~~ ' " ~ 2o^ 



' : exp (~x 2 /2cr 2 



: exp a 



crf(y / a)v27r 

If x ^ [— a\/2a, +erv / 2a] , then 5* 1+ ^ j f ^g^ ) = Vfc, and, consequently, the 
sum obviously vanishes. The fact that Pk{o) is a distribution can be checked by 
integrating over i£l both sides of relation : we obtain 



^Pk (a) = I 

k=0 J-crVZa 



erf (y/a) ctv27t 



exp (~x 2 /2a 2 ) = 1. 



□ 



Consider now a discrete random variable Q (to, a) such that 

(11) Pr (Q ( W , a) = 1 + 1/fc) = Pfe (a) = ^jz^f- ,fc G N . 




12 3 4 5 



FIGURE 1. function / (a) = erf (\fa) 
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As shown on Fig. 1, function / (a) = erf (\fa) gets rapidly close to 1 as a increases 
- typically /(5) = 0.9984 - so that 

, . exp (-a) a k+l / 2 

^ »W = r( fc + 3 /2) 

can be considered as a very good approximation to the distribution of Q(uj,a) for 
large values of a. Thus a stochastic expression of equality (|lt)(l for large values of a 
writes as follows: 

Theorem 1. A Gaussian distribution with variance a 2 can be accurately approxi- 
mated - in the sense of lemma^\- as a q-mixture of Tsallis maximizing distributions 
as defined by 0) where random non-extensivity parameter Q(uj, a) distributed ac- 
cording to 

Figure 2 shows the partial sums X]fc=o£' fc ( a )'-'i+ J - 2a<r 2 i x ) involved in formula 
(|10|l for n = 0, 5, 10 and oo (Gaussian limit), in the case where a = 5 and a = 1. 
We remark that the k — term in the sum corresponds to a uniform distribution 
on interval [— o~\/2a, +er\/2a| 




FIGURE 2. values of J2k=oP k ( a )^i+^,2aa 2 ( x ) for a = 5, cr = 1 
and (from bottom to top) n = 0,5, 10 and + oo 



2.3. The distribution of the non-extensivity parameter. We first remark 
that a classical Poisson distribution with parameter A writes 

(13) q k (X) =exp(-A)^, fceN. 

The distribution p^ (a), as written in H12|l. can thus be considered as a "shifted" 
Poisson distribution. 

Theorem 2. The expectation and variance of a random variable X such that 

Pr{X = k} = p fc (a) 
can be approximated, for large enough a, by 

1 9 
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Proof. The expectation writes 

cxp (—a) fa 
erf(Va) \J n ' \™ 2 



= «?'/, ( " i = — 7 — — \ / - -!- I " : 



^ fcpfc (a 

fc=0 

Thus for a large, m ~ a — 4 . The variance is 

+00 

cr 2 = ^2 k2 Pk i a ) ~ m2 



fc=0 

2 1 ( I'" 

a H a 



4 \ 2 

so that, for large a, a 1 ~ a. □ 

We note the resemblance of this result with the case of Poisson distribution 
for which 

m = A, a 2 = A. 

2.4. Extension to multivariate and power exponential distributions. The 

preceding results can be generalized from univariate Gaussians to multivariate 
power exponential distributions as follows. First we recall the following result from 

m 

Theorem 3. The maximizer of Tsallis entropy with q > 1 under the constraints 

(\Xi\*) = ElXiP* =Ki 

is 

(14) fxAX) = {Mi-T,UW^ 

I else 
where the hyperellipsoid E\ is defined as 

£1 = jxeK"|^A fc |x fe p fc <lj 

and the Lagrange multipliers \k are such that the constraints {1$ are verified: if 
q > 1, all Lagrange multipliers Xk are positive. Moreover, constant A q writes as 
follows: 

rf^-r + i + l") " „ 3,1/7* 



luit/i - = y^? , — . 

Now we are in position to state the following theorem. 
Theorem 4. A generalized exponential distribution 



( n ~, \ X /7fc \ / " 



(Afe > Vfc) can &e approximated as a q— mixture of Tsallis maximizing distributions 
f Xi q (x) as defined by in the following sense: Va > 0, 

B an g (x) Vx S £ a 
e^se 



/Ax,Q(w,a) 0*0 
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where 

A = dtag{al hk },E a = lyeR n \J2x k \y k \^<a 



fc=i 



and Q{oj,a) is a random nonextensivity parameter such that 

Pr (Q (w, a) = 1 + 1/0 = Pi (a) = B a , 7 ex P(- Q ) a ' +T 



is a discrete, "Pois ■son-like" distribution with normalization constant 

-l 



rmalizc 



7 r(i + i/ 7 ); • 

i/erer(i,a) denotes the incomplete Gamma function: 

1 / >+oc i 

r(-,a) = / exp(-t)^ _x tft 

7 ./a 

Proof. Consider random vector y — Ax = la^^x-i, . . . , a i ^ n x n ] T obtained by 
multiplying each component Xi of x by factor a 1 ^ li where a is a positive constant. 
Then 



fvM = ^( 1 --E A ^r ,! ] y y^ E « 

with £„ = {)/€ R n | X)fc=i X kVk k < «} so that 

+°° +°° / ■> !+- A , / 1 " 

piwf A x,i+i(y) = B °,ii,—( — : — 1 --z^ Afc l 
1=0 wr(i + i + i) ^ \ %=i 

( n \!/7fc \ +°° I -a / i n \ 



v 



\k=l v 



□ 

Figure 3 shows the curves that delimit the north-eastern part of the plane (7, 
a) where constant B an belongs to the interval [1 — e, 1] with, from bottom to 
top, e = 10~ 4 , e = 10 -6 and e = 10~ 8 . When (a, 7) belongs to this area, then the 
distribution of the random non-extensivity parameter Q(u>, a) can be approximated 
accurately as 

exp(a a t 

(15) (a) ~ 



r(/ + i + i) 



which is again a "shifted-Poisson" distribution. Its expectation and variance can 
then be approximated as follows. 



Theorem 5. The approximated distribution I15j) has expectation and 

(16) m = a , a 1 = a 

7 



variance 
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Figure 3. areas of the 7, a plane where constant B a belongs to 
[1-e, 1] with, from bottom to top, e = 1(T 4 , e = 1(T 6 and e = 1CT 8 



Proof. Using approximation JTSJl, we get 



+00 



EM0 = («--)(i-(--i)^ I) 



1 r(i-l,o)\ ia^expC-a) 



1=0 

and for the variance 

+00 



r(i) 



a 

7 



1 



r(i-i,«) 



1=0 



all , a 

-3- + ^ + a + a 2 + 2 — 

T 7 7 7 



exp(— a)a^ 1 ( 



1 a 



+ a) — m 



7 7 

a + (--a) 2 ~(--a) 2 

7 7 



□ 



2.5. Corollary: Application to an Hamiltonian system. The result of Theo- 
rem 4 can be used in the following framework: let us consider an Hamiltonian system 
with a phase space of dimension n = 2N. Denote as x = (qi, . . . , qNiPi, ■ ■ ■ ,Pn) 
the phase space coordinates and consider the following Hamiltonian 



H 



k=l 



where are positive constants related to the system. Let us assume that the 
system has maximum Tsallis entropy with parameter q > 1 under the single energy 
constraint (H) = E. Its probability distribution fx writes thus 



(17) 



/x(x) = — l-(«-l)/35>kMZ* 

q \ k=l 



We remark that this distribution coincides with distribution (|14|l , except that con- 
stants Xk in (|14l) are the n Lagrange parameters, whereas there is only one Lagrange 
parameter in l|17|) . namely the inverse temperature (3. However, results of Theorem 
4 apply in the straightforward following way. 
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Theorem 6. The Hamiltonian system with probability distribution JJ7| ) can be 
approximated (in the sense of Theorem 4) as a Boltzmann system with distribution 

g x {x) = Z^exp (^-(q-l)(3j2a k \x k p^j 

with 

and where the non-extensivity parameter q is a discrete random variable following 
the shifted- Poisson distribution f!5p . 



3. Discussion 

We have seen that a Tsallis probability distribution (PD) S q <T i of fixed variance 

(18) S q ^ (x) =A q (l- ^ , 

where the normalization constant is given by JJJ), in conjunction with a shifted 
Poisson distribution p(a) can be used to expand a Gaussian PD (or more generally, 
exponential PDs) of fixed variance in the fashion (for brevity we concentrate here 
upon the Gaussian PD only) 

1 ~ t ~°° 1 / x \ 

(19) —7= exp (~x 2 /2<t 2 ) = Y,Pk (a) —=S 1+1/kA —= , 
cry Ztt cry/ la \cr\/ la J 

Vx £ [—ax/la, +o-\/laj . Now, it is known that for a Gaussian PD the variance a 2 
is proportional to the temperature T (for a purely exponential PD, the temperature 
squared) pQ. Thus, Eq. entails that, for fixed temperature T, a Gaussian PD 
(a purely exponential PD) can be thought of as a superposition of Tsallis' PDs of 
the same temperature. The weights in such an expansion are given by a shifted 
Poisson PD. Comparing Eq. I|19l) with Eq. J2J we detect a striking similarity. Of 
course, in the first we have a "finite" convolution. Moreover, the roles of Gaussian 
(purely exponential) and Tsallis PD's are inverted. 

3.1. Finite thermal baths. The conventional Gibbs' road to derive the canonical 
distribution [2] considers a system S with energy levels denoted by Ci, weakly inter- 
acting with an infinite thermal bath B and assume that one describes the "total" 
system T = S + B by recourse to Gibbs' microcanonical ensemble when its total 
energy E lies in the interval 

(20) E - A < E < E + A, 

with A <C Eo 0. Plastino and Plastino (PP) 12 traverse a new road by assuming 
that B is a finite system, of finite energy Eb and finite number of particles N, 
and in so doing are able to derive, repeating Gibbs' original arguments, not Gibbs' 
canonical distribution but that of Tsallis' (Cf. Eqs. 

(21) k _ Zq . _ Dl _ w _ ? , K|1/<1 -,, 
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Assuming further that, plausibly enough, the bath' number of states per unit energy 
interval grows as a power a — N of E (a common ocurrency [El), PP show that 
the nonextensivity index q* is given in terms of N in the fashion 



Using Eq. © we have then 



as prescribed by the theorems of the preceding Section. 

We can now interpret the contents of Eq. I|19l) in terms of decomposing the Gaussian 
PD in terms of Tsallis' ones. In order for thermal equilibrium to be established 
(and have a Gaussian PD) one requires, of course, thermal contact with an infintc 
bath. This reservoir, in turn, can be thought of as a superposition of finite baths 
containing different numbers N of particles. Thermal equilibrium with any of them 
leads to a Tsallis' PD. 

3.2. Why quasi-Poisson weights in (|19|) ? A possible visualization of our pre- 
vailing state of affairs state is as follows. Due to the short range of the effective 
interactions between our system and the bath (e.g., Lennard- Jones, van der Waals, 
etc.), even if the reservoir is indeed of infinite size the system can not "appreciate" 
such a feature. During a short time interval AT, the system actually interacts 
with a finite number of particles, say M of them. Consequently, we may say that 
the system only "sees" an effective, finte heat bath of size M. Of course, M is 
constantly changing, due to the interaction itself and to the time-evolution. Thus 
our system is constantly "choosing" M particles to interact with out of a uniform 
distribution of them in configuration space. Think, for instance, that the reservoir 
is an ideal gas in a very large box, and that the system, during the interval AT, 
interacts with those particles contained within an effective volume AV. Therefore, 
M would be the number of particles contained within AV. But this is a typical 
situation leading to a Poisson distribution (for the random variable M). Within 
the present scenario we may say that the Gibbs statistics (infinite heat bath) arises 
in an effective way from the superposition (with appropriate weights) of different 
(/-statistics (finite heat baths). In a sense, this situation is the reciprocal of the 
one considered by Beck |3J^|S]. In Beck's formulation, the (/-statistical ensemble 
constitutes an effective description of a system interacting with a heat bath of fluc- 
tuating temperature. Here we show that, under physically sensible assumptions, it 
is possible to regard the Gibbs statistics as describing a system interacting with a 
finite heat bath of fluctuating size. 

It is interesting to examine this mechanism from the point of view of the ergodic 
hypothesis. According to this hypothesis, ensemble averages of dynamical quan- 
tities coincide with suitable time averages 0. Now, if our (/-based origin of the 
Gibbs' statistics is actually realized for a given system, we should expect the er- 
godic condition to hold "differently" in two distinct time scales, (i) In the first place, 
time averages computed during appropriate, finite time intervals AT would coincide 
with (/-averages, with the q- values associated with successive AT-time windows dis- 
tributed in a Poisson-like way. (ii) Secondly, time averages computed during time 
intervals r much longer than AT, would coincide with Gibbs' ensemble averages. 
This kind of behavior has some similarities, as well as important differences, with 
an scenario leading to the (/-statistics that has been suggested by Tsallis' ^U]. Ac- 
cording to this Tsallis' picture, some systems with long-range interactions relax 



(22) 




N -2 
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first to a meta-stable state (characterized by a ^-statistics), which is endowed with 
a long mean-life. But, eventually, the system relaxes to a final state described by 
the Boltzmann-Gibbs' statistics. The mechanism here suggested by us shares with 
that of Tsallis' the fact of (i) being based on two different relaxation time scales, 
and (ii) recovering the Boltzmann-Gibbs statistics in the limit of long relaxation 
times. However, our mechanism differs in one fundamental way from the Tsallis' 
one: in our scenario the system is always described by a g-statistics, if observed 
only during a relatively short interval AT. 

Further developments and applications of our present g-fluctuating approach to 
the Gibbs' ensemble, and its possible verification in (actual or numerical) experi- 
ments, would certainly be greatly welcome. 
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